The goal of this markdown is to determine a conclusion to the genetic architecture (GA) experiment. I have run 4 replicas for each of the GA sets and will summarize the results.
This is a two dimensional co-evolution simulation where newts have developed a toxin and snakes have evolved resistance to this toxin. Resistance and toxicity are the phenotype that I measured and as the snakes and newts co-evolve these phenotypes get larger. As snakes become more resistant, less toxic newts are being devoured, leaving only the most toxic newts. These toxic newts reproduce and create even more toxic newts. The cycle contentiously repeats. However, the growth of the phenotype is not unlimited, there is a burden for having a higher phenotype (fitness cost).
Resistance and toxicity are quantitative traits that are the sum of effects from mutations on these specific phenotypes. In my simulations I create different genetic architectures (GA) for resistance and toxicity to see if there is an effect on the how the population co-evolve. Genetic architectures are a combination of mutation rates and the standard deviation of the mutation effect size. I have ran two types of genetic architecture studies one where one species is at an evolutionary disadvantage (one species will have less variance) and the other study has the same level of variance. In the study where there is no variance disadvantage I test the effects of mutation rate and effect size to see if one or the other has a large effect on how a population can co-evolve. In the study where there is a disadvantage I not only look at the effects of mutation rate and effect size, but also look at how a disadvantage can effect a population undergoing co-evolution. The disadvantage experiment is GA1 (16 simulations) and the same level of variance is GA2 (25 simulations).
Each of the simulations contain their own neutral msprime simulation (to generate genetic diversity) followed by the slim simulation (populations are co-evolving).
GA1 experiment values:
GA2 experiment values:
The first thing I do is gather the data from each GA set, paying attention to all of the parameters especially the replica number. I typically give each simulation a letter associated with the mutation rate and effect size. Here are the trial letters and the parameters associated with it.
| Trial | Snake.Mu.Rate | Snake.effect.size | Newt.Mu.Rate | Newt.effect.size |
|---|---|---|---|---|
| A | 1e-06 | 0.005 | 1e-06 | 0.005 |
| B | 1e-06 | 0.005 | 1e-08 | 0.05 |
| C | 1e-06 | 0.005 | 1e-10 | 0.5 |
| D | 1e-06 | 0.005 | 1e-12 | 5 |
| E | 1e-08 | 0.05 | 1e-06 | 0.005 |
| F | 1e-08 | 0.05 | 1e-08 | 0.05 |
| G | 1e-08 | 0.05 | 1e-10 | 0.5 |
| H | 1e-08 | 0.05 | 1e-12 | 5 |
| I | 1e-10 | 0.5 | 1e-06 | 0.005 |
| J | 1e-10 | 0.5 | 1e-08 | 0.05 |
| K | 1e-10 | 0.5 | 1e-10 | 0.5 |
| L | 1e-10 | 0.5 | 1e-12 | 5 |
| M | 1e-12 | 5 | 1e-06 | 0.005 |
| N | 1e-12 | 5 | 1e-08 | 0.05 |
| O | 1e-12 | 5 | 1e-10 | 0.5 |
| P | 1e-12 | 5 | 1e-12 | 5 |
| Trial | Snake.Mu.Rate | Snake.effect.size | Newt.Mu.Rate | Newt.effect.size |
|---|---|---|---|---|
| A | 1e-08 | 0.05 | 1e-08 | 0.05 |
| B | 1e-08 | 0.05 | 1e-09 | 0.158 |
| C | 1e-08 | 0.05 | 1e-10 | 0.5 |
| D | 1e-08 | 0.05 | 1e-11 | 1.58 |
| E | 1e-08 | 0.05 | 1e-12 | 5 |
| F | 1e-09 | 0.158 | 1e-08 | 0.05 |
| G | 1e-09 | 0.158 | 1e-09 | 0.158 |
| H | 1e-09 | 0.158 | 1e-10 | 0.5 |
| I | 1e-09 | 0.158 | 1e-11 | 1.58 |
| J | 1e-09 | 0.158 | 1e-12 | 5 |
| K | 1e-10 | 0.5 | 1e-08 | 0.05 |
| L | 1e-10 | 0.5 | 1e-09 | 0.158 |
| M | 1e-10 | 0.5 | 1e-10 | 0.5 |
| N | 1e-10 | 0.5 | 1e-11 | 1.58 |
| O | 1e-10 | 0.5 | 1e-12 | 5 |
| P | 1e-11 | 1.58 | 1e-08 | 0.05 |
| Q | 1e-11 | 1.58 | 1e-09 | 0.158 |
| R | 1e-11 | 1.58 | 1e-10 | 0.5 |
| S | 1e-11 | 1.58 | 1e-11 | 1.58 |
| T | 1e-11 | 1.58 | 1e-12 | 5 |
| U | 1e-12 | 5 | 1e-08 | 0.05 |
| V | 1e-12 | 5 | 1e-09 | 0.158 |
| W | 1e-12 | 5 | 1e-10 | 0.5 |
| X | 1e-12 | 5 | 1e-11 | 1.58 |
| Y | 1e-12 | 5 | 1e-12 | 5 |
First, I will look at a summary plot for the different GA sets. There are four points for each GA set (one for each replica).